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ABSTRACT 


This study was an attempt to map the density field in a fin-flat 
plate junction three-dimensionally using holographic interferometry. 
This investigation has extended the density studies by Matulka ie 13] 
and Jagota lied 4] to include, for the first time, interferometric fringe 
information obtained through a transparent model in supersonic flow. The 
fringe information was then inverted by a FORTRAN computer program to 
produce a plot of the density field around the model. The feasibility 
of the method was demonstrated. 

The factors which are thought to have limited the success of the ex- 
periment include vibration of the model, fluctuations of the tunnel flow 
and the fact that the model was somewhat too large in relation to ate 
size of the wind tunnel test section. Schlieren photography was used to 
look through and around the model and to verify that the same flow was 
established as was Be outa by Thomas { 23, 24 | and Winkelman | 26, al 

The data reduction of holographic interferograms was, for the first 
time, accomplished using photographic enlargements. This technique is 
considered to be much easier and more accurate than the one used in the 
previous investigations. However, the data reduction step, because of 
the time and labor involved, is considered to be the rate controlling 


process of the whole analysis. 
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LI. INTRODUCTION 


The determination of the flow field in the wing-body junction of an 
aircraft in supersonic flight presents many problems. The typical ap- 
proach has been to measure the static pressure on the surface using 
many small pressure taps [238 24 | and to determine the velocity field 
in the junction using translatable pressure probes [18, IESE 25). Thomas 
ah, 24 | and Winkelmann { 26, 27 | used azobenzene and oil-smear tests 
to view the flow field around a flat plate-fin junction. From the streak 
patterns on the plate and fin, they were able to illustrate the three- 
dimensional flow field, although in a partially speculative fashion. 
This study has attempted to map the density field in a fin-flat plate 
junction three-dimensionally using holographic interferometry. Although 
this objective was not fully achieved, the feasibility of the method has 
been demonstrated. 

By uSing a Q-switched lees with exposure times of about twenty nano- 
seconds, it was possible to obtain three-dimensional holographic inter- 
ferograms of the density field in the fin-flat plate junction. From 
holograms taken at a number of viewing angles the fringe shifts in dif- 
ferent planes could be obtained. By integrating this information using 
a FORTRAN computer program, the density field can be determined. This 
technique has been previously demonstrated for the flow field of a free 
jet by Matulka (12, 13] and for the supersonic flow field around a cone 
at angle of attack by Jagota [3, 4] , 

The tests were performed at the Naval Postgraduate School, using the 


four-inch supersonic wind tunnel. 





IIT. EXPERIMENTAL APPARATUS 


A. THE WIND TUNNEL 

The investigation was conducted in the Naval Postgraduate School 
blowdown-to~atmosphere supersonic wind tunnel. The test section is 
four inches by four inches in cross section and six inches long with 
two different sets of side walls. The two-inch thick plexiglas side 
walls, which have a refraction index of 1.49, present a complete field 
of view of the flow from the nozzle throat to aft of the test section 
mounting bracket (Figure 1 (a)). The second set of sidewalls used are 
aluminum with high quality optical glass portholes located in the test 
section area (Figure 1 (b)). The interchangeable nozzle for a test 
section Mach number of 2.8 was used for all tests. The nominal run 
time is five minutes at Mach 2.8 with the maximum stagnation pressure 


of about 105 pounds per square inch. 


B. THE HOLOGRAPHIC ARRANGEMENT 

The holographic arrangement is illustrated in Figure 2 and shown in 
photographs included as Figures 3, 4, and 5. The«ajuipment stand was 
rested on a portion of the building floor that was vibrationally isolated. 
A Konrad K-1 pulsed ruby laser with a Pockels cell Q-switching unit was 
used to produce monochromatic light at a wave length of 6943 Angstroms 
and exposure time of twenty nano-seconds. The laser cavity length was 
seventy-three cm. giving a coherence length of about ten cm. To main- 
tain the laser head and output etalon at a constant temperature of 27.5 
degrees centigrade, a Lauda constant temperature circulator Model N was 
used. This was controlled by an electronic relay type R-10 coupled with 


aeedriigan de-ionizer. 





Holograms were obtained by routing the reference beam under the 
wind tunnel and the scemme beam through the test section, and inter- 
secting the two beams on the hologram plate at an intersecting angle 
of approximately 50 degrees. The beam sizes were controlled by trans- 
lating the concave lenses located between the beam splitter and holo- 
gram plate in each beam (Figure 2). The Q-switched laser and optics 
were aligned using a continuous wave helium-neon laser. For reference 
purposes, grids were mounted on the outside tunnel walls and aligned 
using a Ssurveyor's transit. The holographic stand and test section 
were completely enclosed in a wooden box to enable holograms to be 


taken in the daylight (Figure 6). 


C. THE WIND TUNNEL MODELS 

The fin-flat plate models used are shown in Figures 7 (a), (b), (c), 
and (d). The metal portions of both models were aactente cs steel. The 
center section, part of one strake, and all of the other strake of the 
first model in Figure 7 (a) were made of epoxy while the center section 
of the second model, Figures 7 (b) and 7 (c), was fabricated from plexi- 
glas. The flat plate grids in both models were etched into the plastic 
surface and coated with a clear plastic to achieve a smooth surface. 

The models were rotated about their sting mounts as shown in Figure 
7 (d). Alignment for the desired rotation angle was accomplished by 
aligning prescribed lines on the sting mount collar with a scribed mark 


on the sting stand using a surveyor's transit. 





ITIL. ANALYTICAL EVALUATION OF THE DENSITY FIELD 


A. THE BASIC INTERFEROMETRIC EQUATION 

Interferograms are created when two coherent light beams are super- 
imposed on each other and projected on a viewing screen. The light. and 
dark regions observed correspond to the relative phase difference be-. 
tween the two beams which are caused by a difference in the two optical 
path lengths. Consider a coherent beam which is split and then recom- 
bined on a viewing screen. A difference in optical path lengths Pec 
two beams may be achieved in two ways in order to create an interfero- 
gram. The first is to make the physical distance traveled by the two 
beams different. In a vacuum this path length difference is expressed 
as L= (C,4t where ¢C, is the speed of light in a vacuum. The second 
way is to maintain equal physical path lengths but to have the beams 
traverse through different media prior to recombining. In this case 
each light beam will travel at a Bad <2. where n is the index of re-. 


fraction for the medium traversed. The optical path length difference 


then becomes: 
QOL=L (no - ny) = C,At (1) 


The interference pattern or fringes observed may be expressed as a 


function of the optical path length difference or 


4L 


= atl 2 
g > (2) 
where: g = fringe shift 
» = wave length of the light source 
AL = change in optical path 





Combining equations (1) and (2), the fringe shift is then 
g = 4 (ny) - nj) (3) 
» 


The index of refraction is known to be a function of density. Since 
the speed of light is only slightly less in gases than in a vacuum, the 
index of refraction could be closely approximated by the series expan- 


sion (3] 
n= 14+65 (4) 


where 
6 = dimensionless constant related to the Gladstone-Dale 
constant by K = Yes 
(= reference density of 0° C, 760 mm. Hg. 
The variation of with wavelength is small and has a value of 
0.000292 for A = 5893 angstroms. 
Considering a fixed difference in the index of refraction between 


the two beams in Equation (3), then 


s= BE ( GS) ) 


If the density varies in a beam path, the net change in the optical path 
length will be the integrated effect along the beam path or 
L L 
& = &B j CP pees ae) \ f(x,y,Z,)ds (6) 
dfs ray oO 
where 


feo (6a) 
» Ps 


2) 
I 
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(6b) 


Zo.= a plane of constant Z 
ds = incremental distance along the ray 
In order to determine the density along the beam path where the fringe 


shift is known from an interferogram, Equation (6) must be inverted. 


B. THE INTEGRAL INVERSION 

The integral inversion technique was ee reported by C. D. Mal- 
donado et al in 1965 cr 10, at R.D. Matulka (blag 13| ang sh sae 
Jagota [3, a] used this method to determine the density variation in 
an asymmetric free jet and about a cone at angle of attack, respectively. 
The technique involves representing the function, f(x,y, %) in Equation 
(6) by a complete set of orthogonal functions where the unknown coef- 
ficients are evaluated using the orthogonality relationship between the 
set of functions. The functions are orthogonal over the entire plane 
and also have the propery of being invariant in form to any rotation of 
the coordinate system. Figure 8 illustrates the coordinate system for 
the inversion where x and y are the fixed laboratory coordinates and x' 
and y' are the coordinates in which the fringe number function is defined. 
As the view through the test section is varicd the primed coordinates 
are rotated with respect to the fixed coordinates x and y. 


The fringe shift expressed in Equation (6) may be written as the 


transform 


BO SRY 3) 2, = ERK Ze) (7) 


I: 





or, inverting the equation, the density function, f£, is equal to: 


~{ 
E(x, ¥, 20) = | e(&. ¥ > 2) (8) 
The density function can be expanded in the following manner using a 
=m 
set of polynomial functions, U a (4X, AY), and unknown comples 
+ 


cw 
Becetmenents, Ci woe (A) 


FF cal CN MOU leer? Cra Cae] 


Oo Kk-0 


where a ¥2, ™M=O 
1 = mel,2,5,""° 


A= arbitrary scale factor 


£(x,y¥5z,) = (9) 


ve 4 *) 
Cc 


The polynomial functions, Ua @xX,Ay), are invariant in form to a 
rotation of the coordinate system [10, LZ Us|): They also have a Gauss 
transform which makes them adaptable to the physical situation and to 

manipulating into the form of Equation (7). The functions are defined as; 


=m k OP KI(AK eATYZI™ 1 So times pm 
X,Au)= (- : pee z a 
Sta , 2) ( ) «| Time K)! | Sc Is (« K+ d*y* ) (10) 


where b = tam ( ¥) 6 i (10a) 


Es = Laguerre polynomial 


k 


; a f nae 3)) omen (ax trary) | ae 


And the Gauss Transform of US, LS: 
£m a ae, sims (Ay') (41) 
' - - & _Hine2k sm 
if (ay58)=| U™ (ax aye” dx's aero 
me 2k ioe meZk [! (m+k)! |* va 
where (xu')= Hermite polynomials 
Hine zk 3 


By applying the transform above to Equation (9), the fringe function in 
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Equation (7) can be written as: 


£2 mm ~Eme -a'y'2 
a,( $y 412c) <) \ : ne ae” ‘ PTE Core J How en’e (12) 
{ kL msk)! Dare 


using the following orthogonality relationship on Equation (12) 


tT tim€ ting hry ; 
\ A A as | [4 awakes Anes (tye J dy = 





mero Kro 
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¥, 
2K ig 
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where 4 is the kroneker delta, the expansion coefficients C mews 


can be determined by: 


‘ 7) 
tm ox [(ki(msk)!)2 ; Ly 
A) 2—, |—- , (14) 
Cea! } S| (m+ Zk)! 9(4, 5,26) Se (eye. i a§ 
= 
Substitution of the coefficients in Equation (14) bact into Equation 
(9) results in the density variation being expressed as: 


co © y ~ (A*K™ 4A*u*) 
L(x, eee )uc (ay) ) ém[ kL Omani] = ‘ 


mzO K20 


ff coo i 

REAL | | | gtuisz0e ms Ames 4")du 5 | Doe (AX,h uy) (15) 
- =a 

or by inserting Equation (10): 


K mm m 
(-1) KI] yboueanes VA |. (h*x%4* 2). 
$x,y22= (F ays ) & ™“(ms+zk)! ( ens és - J 


mx0O 


(Amy? Au? 
[ BT (4) coslm$) * Dee (4) suntmg)] xr ear y*) (16) 


where: 


Pe (a) f I ay',§, 2.) Cos(mg) H Pier yd € (17) 
a 


Dirvak (4) f' [ gcyiseoou (m€) H oe Au’) dy d§ (18) 
= Tf i 


is 
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Equations (16), (17), and (18) are the basic equations used to 
calculate the density distribution from the experimentally determined 


fringe variations. 


C. THE NUMERICAL PROCEDURE 

The form of the density distribution in Equations (16), (17), and 
(18) must be modified in order to input the experimentally determined 
fringe distribution. First from Figure 8 and Equation (6b) it can be 
seen that it is only necessary to integrate Equations (17) and (18) over 
an area where the density is changing from a known density, Peo . Out- 
side of this region where there is no change in density, the function 
E(x,y4, - 0, 1.€. outside the test section. Also since the fringe 
distribution is taken in small increments over the test area the coef- 
ficients, B and D, can be approximated as 


oe Sua Xie 
Bm ws). Fa c. +A‘, sy %, + DX;) | ec Cos wnsrds | Wome ag *rrda (19) 
J 


é<l yso 


tr! j- S 
Di “) oer AK; ) sini | Ey seen) dx (29 


C= jeo xi 
The integral of §$ is easily aeeeamucd and by using the derivative for- 


mula for Hermite polynomials the integral of x may be manipulated to yield 


=! 
a “ote ime awn » S als, +h%, s%i + QAXx;) ¢ 


£20 sso 
Lsun(m§,,,) som m8) 1LHaaasist rH aone ih XO (21) 
Disk (A) © - Laxmanetnn | 3 PG 105, xi + Oxi) © (22) 


[eos (m$, ,)- cos€m§. STH, ons (AX) ~ Hine oko AX) J 


In the Pe anaievion of the density function from Equation (16), 0 


taining the infinite summations experimentally is not plausible or 


14 





possible. It has been demonstrated that by using a finite number of 
terms and by adjusting the values of AS , & x and A, it is possible 
to obtain the density distribution with very good accuracy oe ae a 


13]. Equation (16) then becomes: 


K 
= ay) Seale i)" ki == *u*) e 
boxy,22 (=) £— L_. | (A*x*+ A*y*) 


(23) 
" (Arx*+A'y*) [ ey (A) COs(ms$) + Da) sun (eng) | m~ 
~(K*X* + KU) 


IV. EXPERIMENTAL PROCEDURE 


A. LABORATORY TECHNIQUES 

The analysis of a free jet by Matulka [12, 13] illustrated how 
holographic interferometry can be used to obtain a complete three dimen- 
sional plot of the density within a moving transparent flow field. 
Jagota ee 4 | in his study of a cone at angle of attack in supersonic 
flow went one step further by introducing an opaque object into an 
assumed steady state flow field and describing the density field three- 
dimensionally. 

This investigation has attempted to determine the three-dimensional 
density field around a transparent object in a supersonic flow field by 
pasSing a light beam through the object. Specifically the interest was 
to describe the flow field existing in the junction of a fin-root inter- 
section. 

1. Model Considerations 
In order to obtain uniform flow around the fin-root area, a model 


of the form shown in Figure 9 was selected. The flat plate has a knife 
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edge and is intended to remain at zero degrees angle of attack so as 

to establish the flow conditions illustrated. The fin edges were made 
circular in order to approach flow conditions similar to those estab- 
lished in Winkelmann's [ 26, 27 | investigation. Plastic and metal strakes 
were added on the model sides so as to maintain two-dimensional flow as 
well as to add strength to the flat plate. In the first model constructed 
(Figure 7(a)), maximum visibility of the plastic fin-flat plate center 
section was achieved by bolting the leading edge and aft plate together 
through the oleae center section. Due to model flexure, this design 

was found to be unsuitable and the second model in Figures 7(b) and 7(c) 
was constructed. The model was made from a single piece of stainless 
steel. The model rigidity was satisfactory but unfortunately the strength- 
ening borders around the plastic center section reduced the holographic 
visibility somewhat. 

Since the wind tunnel blocks were fixed, it was possible to make 
multiple test runs with the same flow conditions over the model provided 
supersonic flow had been established over the model. 

2. Holographic Techniques 

In order to obtain holographic interferograms it was necessary to 
ensure that the optical path lengths of the scene and reference beams 
remained approximately equal. Since the ruby laser is believed to have 
a coherency length of approximately ten centimeters, the quality of 
lengths is far less critical than in the classical Mach-Zehner inter- 
ferometric approach. Consequently a string was used in the experiment 
to trace the reference beam and then adjust the scene beam. This method 
kept the two beam lengths within one centimeter of each other. Since the 


scene beam traversed approximately 4.5 inches of plastic tunnel walls and 
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grids whieh the referenee beam did not, it was neeessary to eompensate 
by making the scene beam physically 2.25 inches shorter than the re- 
ferenee beam. The referenee beam varied in length from 61 inehes to 
68 inehes during the experimentation. 

To determine the fringe/density field, finite fringe interfero- 
grams were made by three different techniques. [In the direet fin-root 
flow approaeh the diffuser plate was part of the model. Either the 
Mirror , Ms » in Figure 2 or the hologram plate’ holder was translated be- 
tween the no-flow exposure and the flow-established exposure. [In the 
total model flow method the diffuser plate was loeated between seene 

“ beam lens, L3, and the test seetion (See Figure 2) and it was translated 
horizontally or vertieally. Translations were varied from .001 inehes 
to .006 inehes with the translation distanee of .003 inehes yielding the 
best fringe separation. | 

Most of the holograms taken using basically the holographie ar- 

rangement shown in Figure 2 gave well-defined fringe patterns. In order 
to improve upon the fringe definition, a variety of teehniques was at- 
tempted. The transverse mode selector was varied from 1.0 mm. to 2.5 mm. 
in inerements of 0.5 mm. to determine the best lighting of the model. 
The hologram plate holder was rotated horizontally to various positions. 
These positions varied between being perpendieular to the seene beam to 
being perpendicular to the bisection of the angle between the scene and 
referenee beams. Polarization plates were added in both the seene and 
reference beams between the test seetion and hologram plate in the seene 
beam and between the last mirror, M4, and the hologram plate in the re- 
ference beam. A one-quarter wave plate was also placed between the first 


lens, L,, and the beam splitter as reeommended by Okayama and Emorl (14]. 


17 





The holograms were taken using 4" x 5" Agfa-Gavaert 8E-75 holo- 
gram plates. The developing process involved: 

1. Five minutes in Kodak D-19 develcper 

2. Thirty seconds in an acetic acid stop bath 

3. Five minutes in standard fixer 

4. Five minutes in a flowing water bath 

5. One minute emersion in Kodak Photo Flo wetting agent 

6. Drying using blowing cool air 

a. Direct Fin-Root Flow Method 

Since the flow field in the fin-root junction is assumed to 

be identical on either side of the fin, then it is only necessary to 
determine the density on one side of the fin. To accomplish this it is 
necessary to obtain holograms for 180° of view as shown in Figure 10. 
The holograms for the views from 0° to 90° can be obtained by using a 
frosted fin and flat plate as shown in Figure 11. The advantage of 
having the frosted plate as part of the model is that the fringe/density 
information obtained by the interferogram is believed to be only for the 
area between the fin-root intersection to the tunnel wall vice the whole 
test section, but this was not verified. In order to obtain the fringe 
information for angles greater than 90° but less than 180° the fin 
would be exchanged for one containing a stainless steel reflective sur- 
face. The scene beam would then enter the test area from the viewing 
port below the tunnel and be reflected to the hologram plate as shown 
in Figure 12. 

b. Total Model Flow Method 

In this method the diffuser plate was located in the scene 


beam outside the test section as shown in Figure 14 and the flat plate 


18 





center section and fin were made of optically clear plexiglas. By 
translating the diffuser plate between exposures of the hologram, an 
interferogram of the whole density field in the test section about the 
model was obtained. From Figure 13 it can be seen that due to symmetry 
only 90° of view was required to obtain the density field. This makes 
it much easier experimentally to take the holograms than the previous 
method described. 
3. Schlieren Analysis 

A standard Schlieren knife-edge system was used to verify the 
establishment of the supersonic flow network around the model as shown 
in Figure 9. Photographs were also taken of the flow with the model at 
0° and 90° rotation in order to compare the fin shock conditions with 


those obtained by Winkelmann [26, 27]. 


B. PHOTOGRAPHIC TECHNIQUES 

By illuminating the Eoloetane ale a helium-neon laser beam which 
has a wave length of 6328 Angstroms, the original scene was reconstructed. 
Since the original scene beam and the reconstructed beam were of different 
wave lengths, there is actually a small distortion in the reconstructed 
scene but of neglectable effect because the hologram plate emulsion in 
the development process also shrinks. 

The typical method for reconstructing the scene is to il lumimatesthe 
hologram as illustrated in Figure 15. The qarnuee glass used in the con- 
struction of the hologram appears to act as an infinite light source of 
non-parallel rays which illuminates the scene. If a small aperture 1s 
positioned at the focal plane of the imaging lense, an almost parallet 


set of rays may then be selected as shown in Figure 16. A third method 
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illustrated in Figure 17 uses a small diameter conjugate beam to illu- 
minate the hologram. A large depth of field is achieved because the 
narrow beam acts as an aperture. This effect was of considerable ad- 
vantage since it enabled both the front and rear grids, the model, and 
the fringe patterns to be simultaneously projected on the screen. The 


best photographs were obtained by focusing on the plane of the fringes. 


C. DATA REDUCTION 

Photographic .interferograms were obtained by illuminating the holo- 
gram scene with a thin laser beam and using a camera with a viewing screen 
located in the film plane as shown in Figure 15. The line of sight in 
the plane desired was achieved by translating the hologram until common 
points on the front and rear grids were aligned. The camera, with the 
aperture set wide open at 7.7, was then adjusted to give the best focus 
on the fringe plane. The best photographic results were achieved by 
using an exposure time of 1/10 second with Polaroid Type 55 P/N film. 

It was felt that the density field could be well defined three-di- 
mensionally along the fin if the density fields were determined in four 
planes perpendicular to the Z-axis andqjually spaced along the fin as 
shown in Figure 18. [In obtaining the fringe data across a constant Z- 
plane it would be necessary to take six photographs per rotation angle, 
aligned at appropriate intervals down the y' axis, and then graphically 
mate the fringe data to form one complete set. The six photos across the 
field were felt necessary because the optical path length from the model 
to the hologram plate varied for those points not on the aligned plane. 

The fringe shift reduction was accomplished using two different tech- 


niques. The first was to project the negative, using a photo enlarger, 
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onto a sheet of paper and trace out the fringe pattern, model surfaces 
and grid lines. The light fringes were traced out since it was much 
easier to judge their center line. From the fringe lines forward of 

the fin in the region of uniform flow one fringe line which appeared 

the straightest and compatible with most others was selected as the 
bench mark. A straight line was then drawn over that fringe and extended 


past the aligned Z plane (i.e. y axis). Lines parallel to the bench 
mark line were then drawn likewise over the remaining fringe lines. The 
fringe displacements were then read relative to the lines drawn at the 
points of intersection of the fringes with the y' axis. The radius of 
the inversion circle was selected so that the fin-root intersection was 
the origin and the fin tip was the 100 percent point. 

In the second technique an enlarged positive photograph was made 
from the negative. Again one fringe line in the uniform flow region 
just forward of the fin which appeared to be parallel with the majority 
of the fringes was selected as the bench mark. The remaining fringes 
were likewise traced over with lines parallel to the first. The fringe 
displacements were then measured relative to the lines drawn. For fur- 
ther details see Appendix A. 

The locations at which reference lines crossed the y' axis in both 
techniques above were further adjusted to account for the tunnel wall 
refraction displacement as shown in Figure 19 and computed by a computer 
program in Appendix B. Once these corrections were made, the radial 
variation of the fringe number could then be plotted for the various y' 
alignment planes and a smooth curve drawn through the data points. The 
fringe number at 201 equidistant points across the field can then be 


obtained for input into the computer program, HOLOFER, in MODE 3. For 
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further details on how to use the computer program, HOLOFER, see Appen- 
dix C. Once the data from all the rotation angles of the model have 
been put into the computer program, the program will then calculate the 
density field across the inversion circle for that Z-plane. After the 
density has been calculated for all four Z-planes in Figure 18, a three- 
dimensional plot of the density field can be made by connecting points 


of equal density across the fin. 


V. EXPERIMENTAL RESULTS AND DISCUSS ION 


The initial attempts to establish uniform flow over the flat plate 
shown in Figure 7(a) were unsuccessful due to model vibration and flexure. 
Movement of the model sting within its holder and flexure of the model 
plastic center section allowed the model leading edge to establish a lit- 
tle over 1° angle of attack upward when flow was established. Due to 
various modifications made in attempts to eliminate the vibration, the 
sting finally fractured. 

In an attempt to eliminate these problems, the second model shown in 
Figures 7(b) and 7(c) was made of a single piece of stainless steel and 
the sting was mated to its holder to within .001 inches. The model center 
section was made of poured epoxy and the strakes were both made of stain- 
less steel with plexiglas inserts. The model rotation about the sting 
Was reduced to approximately 0.3° angle of attack upward. Due to the 
vibration in passing through the transonic range and a weak glue seal 
between the metal strakes and plexiglas inserts, the inserts were found 
to break loose. They were subsequently removed and not replaced. 

In order to determine the flow field using the direct fin-root approach, 


the epoxy fin and model center section were frosted on one side using fine 
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emery paper (see Figure 20). Holograms were taken of the model at ro- 
tation angles of 0°, 12°, 45° and 90° using the holographic arrangement 
in Figure 2 excluding the diffusery plate between lens, La, and. the test 
section. The mirror, Ms » Was translated in various directions from to- 
wards to parallel to the test section between the aaa ae flow 
and with flow established. Fringe patterns were obtained around the 
model, but only at 0° and 12° rotation could any fringe patterns be ob-- 
served across the fin. It was found that the fin fringe pattern appeared 
to remain almost unchanged no matter how or how much the mirror, Ms, Was 
translated while the fringe around the model changed appropriately... For 
instance, in Figure 21 the mirror was not translated and in Figure 22 the 
mirror was translated .006 inches horizontally parallel to the tunnel. 

Since fringes could not be observed across the flat plate at 45° and 
90, it was felt that the epoxy center section might be too imperfect 
optically. Consequently double-exposed holograms were taken with no 
flow through the test section and various translation distances from 0 
to .005 inches. The fringe patterns were excellent across the whole 
model and their spacing decreased according to the increase in the trans- 
Pation of Ms. Next the diffuser plate in the scene beam in Figure 2 was 
inserted with the model at 0° rotation angle and translated between no 
flow and flow exposures of the hologram. The fringes about the model 
were of excellent quality but the double diffusion of the scene beam 
through the model caused all fringe patterns on the model (Ein) pico 
disappear. 

It was felt at this point that the fringe pattern obtained across the 
fin was caused by the movement of the model to an angle of attack and 


possibly by model vibration, although none was observed visually. The 
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lack of any fringes across the flat plate center section is not well 
understood but is believed to be caused by vibration of the model. 

Since it was not possible to obtain acceptable interferograms with 
the diffuser plate as part of the model due to model motion, it. was felt 
that the effect of minor model movements could be eliminated by using 
an optically clear model and an external translating diffuser: plate. 
Therefore the fin and model center sections were replaced with optically 
clear plexiglas. With these changes it was found that the flat-plate 
leading-edge ee of attack had been reduced to approximately 0.1°. 

Initially double-exposure holograms were taken using the arrangement 
in Figure 2 and translating the diffuser horizontally. At O° model ro- 
tation the holographic interferograms were excellent. But as observed 
im Figure 23 it would be extremely difficult to determine the fringe 
change across the fin since no free stream reference fringes were avail- 
able, due to the flat plate leading edge Prandtl-Meyer expansion and fin 
shock intersecting the tunnel top just above the fin. With a larger tun- 
nel or smaller model this would be an excellent technique. 

The diffuser plate was then translated vertically between exposures 
and excellent horizontal fringe patterns were obtained as shown in Figure 
24. The model was then rotated to 223° and the same holographic technique 
was used. Excellent fringe patterns were obtained above and below the 
flat plate center section. Fringe patterns across the flat plate center 
section and fin in this area were very light and usable interferograms 
could not be photographed with the polaroid camera. Im an attempt Co 
improve on the fringe quality, polarizer plates were inserted in the 
reference beam between the mirror, My and the hologram plate and in the 


scene beam between the lens, L., and the diffuser plate in order to ensure 
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polarization of both beams. No significant improvement could be noticed 
and they were subsequently removed. A one-quarter wave polarizer plate 
was then placed between the lens, L, and the beam splitter in order to 
utilize circularly polarized light for the reference and scene beams. 
Okayama and Emori [14] found that their image resolution improved con- 
siderably; however, with this particular arrangement little to no improve- 
ment in the fringe resolution was observed and the approach was abandoned. 

The transverse mode selector was then varied from 2.0 mm. to 1.5 mm. 
and later to 1.0 mm. in an effort to improve the coherency length of the 
laser light and consequently the image resolution. Due to the decrease 
in output light intensity, up to six exposures were taken during a run 
with flow established. Image and fringe definition were not found to 
increase possibly due to model vibration which was not visible to the eye. 

The model was rotated to 45° and 673° and double exposure holograms 
were taken with a 2.5 mm. transverse mode selector. There were no ob- 
servable fringe patterns in any portion of the test section. Consequently 
it was believed that supersonic flow was not established due to tunnel 
blockage caused by the shock wave from the model and by slight model 
vibrations in passing through the transonic range. In the transition to 
supersonic flow, the model leading edge would sometimes flex as much as 
a2 depending upon the transition time. 

The test section walls were changed from the total plexiglas side walls 
shown in Figure ea) to the aluminum walls with the optical quality glass 
port holes (Figure 1(b)). The better quality glass would hopefully im- 
prove the viewing and the port holes made the model much more accessible. 
The metal strakes were also removed from the model since they appeared 


to have a minimal effect on the flow, were an interference optically, 
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and appeared to add very little structurally. Flow tests were run and 
the model was found to attain an angle of attack of about 0.14° upward 
when flow was established at Mach 2.84. The angle of attack was tol- 
erable, however the model also maintained a visible high frequency 
vibration of about + 0.05°. Double exposure holograms were taken with 
the model at 0° rotation and fringes could only be observed Ota NOE 
the Prandtl-Meyer expansion on the top of the model and forward of the 
flat plate leading edge shock on the bottom. The tunnel seals were 
found to leak and were replaced and resealed. The same angle of attack 
and high frequency model vibration still persisted. Double exposure 
holograms were again taken and a very faint fringe pattern was observed 
across the model. In reproducing the scene when aligned along a constant 
Z-plane passing through the fin, it was impossible to photograph the 
fringe pattern across the fin because of a bad speckle effect. Double 
exposure holograms were made at model rotation angles of 45° and 90° with- 
out and with flow. The holograms taken without flow gave excellent fringe 
patterns across the model however no fringes could be observed with flow. 

At this point a knife edge Schlieren system was set up to determine 
just what the flow conditions were in the tunnel and around the model. 
The first runs were made at a model rotation angle of 90° and presented 
the shock knee in front of the fin (Figures 25(a) and 25 (b)) which was” 
pointed out by Thomas | By 24 | and Winkelman | 26, 27]. Due to the model 
vibration, the shock extending out from the shock knee was observed to 
flutter about +% cm. 

The model was then rotated back to 0° and Schlieren ciocomeene were 
taken as shown in Figure 26. During all but one of the runs at this 
rotation angle, the model would pitch up approximately 0.14° angle of 


attack and continue to vibrate about + 0.05° at a high frequency when 
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flow was established. The plate leading edge shock and fin shock 

appeared quite fuzzy and Pent aeiuc to O1f efron the tunnel Supply reser- 
voir mixing in the FLOWS a Itene oi] smear pattern can be observed across 
the fin in Figure 26. During one run, which could not be duplicated, the 
model pitched up as flow was established but did not vibrate. The plate 
shock and Prandtl -Meyer waves, fin shock system, and tunnel Mach lines 
became very distinct and well defined as seen in Figure 27. The schematic 
of the Schlieren photographs in Figure 28 nie out the cause for the 
various flow lines observed. The non-uniformity of the free stream caused 
by tunnel leakage can be easily seen. Consequently the experimental work 
was discontinued due to tunnel conditions and time considerations. 

It was felt at this point that if the holographic interferogram taken 
of the clear plexiglas fin at 0° rotation angle could be reduced to use- 
ful data for the computer program then the holographic method would be 
Co a certain extent verified even though the actual density field could 
not yet be determined. 

In obtaining interferograms from the hologram, the reconstruction 
technique shown in Figure 15 was used. The plane of constant 7 across 
the model to be reduced was chosen to be the forward most Vegtical grid 
line crossing the fin. It should be pointed out that in order to simplify 
the hologram alignment process, all four reduction planes across the model 
Should have been scribed on the exterior grids. The three photographic 
pOints across the Z plane on the y' axis were, for convenience, chosen 
to be where horizontal exterior grid lines crossed the y' axis on the fan. 
Five interferogram photographs were taken at points A, B, and C as shown 
in Figure 29. Photographs 2 and 5 were taken uSing a Kodak Wratten Gelatin 


Filter N.D. 2.0 placed in the reconstruction beam in order to reduce the 
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beam intensity and increase the fringe definition. Thcse two photo- 
graphs were also used to provide a check on the consistency of the re- 
duction process. 

All five negatives were blown up on the photo enlarger in the dark 
room and the plate, fin, grid lines and fringes were traced out on a 
sheet of white paper as shown in Figures 30-34. It was very difficult 
EeGmtracemthe fringes in the region of the fim ta@p and fin root due to 
the photographic resolution and the fin tip shadow. Also fringes were 
not visible in the region of the fin leading edge shock. Consequently, 
connecting the fringe in front of the fin to the correct fringe on the 
fin was a best guess effort. Typical of the problem was that fringe 
lines in two of the drawings were initially improperly connected across 
the fin leading edge shock. After bcing checked against the photograph 
negatives,the fringes were reconnected and the data taken correlated 
well with the data from other drawings. The fringes in the fin root 
area were extremely light in some photographs which madc it quite dif- 
ficult to determine their centerline crossing the y' axis. Another 
point of difficulty was determing the location of the top and bottom 
of the fin. An error in drawing here will have an effect later when 
the fringe locations are normalized with respect to the fin height. 

Enlarged positives were then made from the negatives as shown in 
Figures 35-39 to see if the data accuracy could be improved by eliminating 
the difficulties of tracing the fringes in the dark room. This method 
also made it possible usually to reassess the assumed path of the fringe 
lines at a later time if the data did no appear to correlate properly. 
The same problems of locating the fringe lines in the region of the fin 
tip and of connecting the fringe lines across the fin shock are readily 


apparent in the figures. 
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To obtain the fringe change across the fin one of the straightest 
fringes in the free stream region which paralleled the majority of 
other free stream fringes was selected. A straight line was drawn 
through its centerline and extended to cross the y' axis. The other 
free stream reference lines were then drawn parallel to the first and 
adjusted so as to best follow the centerline of the selected fringe. 
This method actually averages the free stream conditions and is only 
valid if the free stream is essentially uniform. Sinceit was not pos- 
sible to obtain in the photographs the free stream fringe pattern for- 
ward of the leading edge Prandtl-Meyer expansion for the lower fin region, 
the fringe reference lines in photo enlargements were drawn over the 
fringes just prior to the fin shock. In order to adjust them to free 
stream conditions their location was moved downward a distance equal to 
1.1 times the average fringe interval for that photograph. The figure 1.1] 
was observed in all the enlarged photographs to be the approximate fringe 
change across the Prandtl-Meyer expansion. In the drawings (Figures 30-34), 
fringe reference lines for the lower fringes (generally y'4.8) were 
initially aligned along the straightest portion of the fringe prior to 
intercepting the fin shock. By comparing the reference line location on 
the drawing with the fringe pattern in the photographic enlargements the 
reference line location and fringe change were adjusted by an appropriate 
portion or all of the 1.1 fringe change caused by the Prandtl-Meyer ex- 
pansion. All of the reference lines were then corrected for the tunnel 
wall and grid plastic parrallax shown in Figure 19 and computed in Appen- 
dix B. The reference fringe locations were then normalized with respect 
to the wing height and the fringe numbers calculated by dividing the 
fringe change by the average fringe interval. For further details and 


calculations see Appendix A. 
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The data obtained for alignment points A and C by the two different 
reduction processes are plotted in Figures 40-43. The data obtained 
from the photo enlargements aligned at point C and shown in Figure 43 
gave the best data agreement between two photographs. The worst data 
agreement was obtained from the reduction of the drawings aligned at the 
same point. The inconsistency in the data was probably caused by con- 
necting the wrong fringe lines across the fin leading edge shock. The 
data fluctuations and discrepencies between the curves in the figures 
could have been caused by a number of things. [It could have been caused, 
for instance, by not drawing the fringe reference line parallel to or 
exactly on the free stream fringe center line or slightly missing the 
fringe center as it crosses the y' axis or by misconnecting fringe lines 
across the fin shock as was illustrated in Figure 42. Fringe location 
errors could be caused by misdrawing the fin tip and root lines as was 
mentioned earlier. Another contributer would be measurement errors. 

To analyze these sources for error, first consider that the typical 
fin size in the drawings and photo enlargements averaged about 2.5 inches 
high and 2.75 inches wide and that the fringe spacing averaged about 
0.10 inches. All measurements were taken using a ruler graduated in 
0.01 inches and readings were made to the nearest .005 inches. Since the 
average difference between the data points and curves ran around ¥% fringe, 
some figures were calculated to determine what measurement errors could 
produce this fringe error. It was found that an error of :025 in¢ehes in 
alignment of the fringe reference line with the free stream fringe center 
line and/or fringe center line crossing the y' axis could produce ¥ fringe 
error. This fringe error will also occur if the fringe reference line 
differs from the free stream center line by more than 1.4° when drawn 1 


inch from the y' axis or by more than 0.36° when drawn 4 inches from the 
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y' axis. A difference of .01 inches in the average fringe interval 


could also produce a %-fringe error; however, this is not too likely 
Since it is an average of fifteen to twenty-five intervals. By com- 
paring the actual height-to-width ratio with those found in all the 
drawings and photo enlargements it was found that average error was 
around 2% or a distance of .02 on the y'-axis. 

In order to compare the interferogram negative quality and the two 
different reduction techniques, a plot of the data for each negative was 
made as shown in Figures 44-48. Photographs 1 and 3 in Figures 44 and 
46 gave the smoothest curves indicating the highest interferogram re- 
solution. Photograph 5 (Figure 48) gave the worst dispersion indicating 
poor interferogram resolution; yet looking at the enlargement in Figure 
39, the fringe line contrast is very good. In comparing Figures 40-48 
it appears that the reduction technique using the photographic enlarge- 
ments gave the most consistent data. This technique also provides a 
much easier and faster recheck on the proper tracing of fringe lines 
because it is extremely difficult and tedious to duplicate. the fringe 
pattern to the same scale over the drawings using the photo-enlarger in 
the dark room. 

All the data obtained by either method for one alignment point were 
then plotted in Figures 49 and 50 in order to observe the dispersion. 

A fringe number dispersion of about 0.6 fringes was observed in the data 
taken from the negatives aligned at point A and a dispersion of about 
0.8 fringes for the data about point C. The dispersion is attributable 
to the inaccuracies in the drawings, to the photographic quality of the 
interferograms, and to the inaccuracies in correcting the lower fringe 
reference lines to free stream conditions. The last point is based upon 


the increased dispersion between the fin tip and fin root data. 
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Figures 51 and 52 show an integrated curve of the fringe change across 
the wing as determined by cach reduction method. In constructing the 
curve, the data from the three aligned points was plotted so as to just 
overlap each other. These two curves were then compared in Figure 53. 
The maximum variation in the fringe number is about 4 a fringe but the 
Variations in the location of the fringe maximums and minimums average 
about Q.1 inches on the fin. The location difference is probably due to 
improper drawing of the fringe reference lines compounded with not being 
able to measure the wing height accurately. Due to the inaccuracies 
introduced in tracing the negative and then reducing the data it is felt 
that the photo enlargement method is the more accurate reduction method. 

With fringe data from only one field of view, the density field, ob- 
viously, could not be obtained. It was felt useful to consider the 
flow field axisymmetric in order to exercise the computer program and to 
provide a check on the program's ability to handle these particular curve 
shapes. Fringe data at 101 equidistant points across the fin from 
OS Y'S 1.0 was obtained from Figure 53 for both curves and fed in HOLOFER 
in Mode 3 for the axisymmetric case. For further details on HOLOFER sce 
Appendix C. The scale factor, & , in Equation 9 was then varied from 0.2 
to 2.5 and a value to 1.0 was determined to ea the most accurate den- 
sity solutions. This value was verified by feeding the function data, 

( eco -1), calculated by Mode 3 back into the program in Mode 1 and 
comparing the fringe data calculated with the original fringe data obtained 
ipom Figure 53. 

The density distribution for both the ‘ems and photographic reduc- 

tion cases is plotted in Figure 54. The density as Y' approaches zero 


actually goes as the drawn lines even though the points imdicate’ a%dip. 
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Matulka [ 12, 13] pointed out that the computer program accuracy does 

not converge at the origin. The large variations of the density curves 

at values of Y') 0.8 are caused by the program trying to adapt to a step 
or shock wave type function at y' = 1.0. If the remainder of the fringe 
data in Figure 53 for values of Y' > 1.0 had been included as input data 
the density curves would have smoothed out. This shock wave step function 
effect was demonstrated and analyzed by Matulka f12, 13). 

The low values of density for the photographic data curve around 
Y' = 0.38 were unexpected but not surprising. First, the photographic 
fringe data curve is more extreme than the drawing fringe data curve and 
second, the density curves are not true values anyway since the field 
was considered axisymmetric and this is not the case in reality. 

In general the density curves are felt to be reasonable under the 
assumed conditions. Consequently it is believed that the computer pro- 
gram could very easily and accurately handle a complete analysis of the 
flow field around the fin-flat plate. 

In completing the analysis there are some data reduction problems 
Which would have been encountered in reducing the interferograms at 
other model rotation angles which merit discussion. With no model ro- 
tation, adjusting the fringe reference lines close to the fin root to 
account for the leading edge Prandtl-Meyer expansion was relatively easy. 
However, when the model is rotated, the fringe shift across the Prandtl- 
Meyer expansion can no longer be considered a constant and it will be 
very difficult to adjust the fringe lines at lower values of Y' to free 
stream conditions. The best solution would be to increase the scene 
beam diameter and/or reduce the model size in order to photograph the 


free stream fringe lines forward of the plate leading edge. The free 
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stream lines could be connected to the appropriate fringe lines forward 
of the fin and this would eliminate the numerical correction and increase 
the fringe data accuracy. 

Also at the rotated angles, fringe information will not be available 
for portions of the reduction plane due to shadows cast by the. model 
Sides as shown in Figure 55. In the lower portion of the hologram the 
fringe curve can be connected with a smooth line because the density 
field in that portion should be essentially constant. Care must be taken 
in completing the curve, though, since this information will be used in 
the integration of other rotation angles. The fringe curve in the upper 
portion must also be completed carefully tut should be easier since the 
shadow will be smaller. For this particular model it was calculated 
that for angles greater than 46.6° the upper shadow would not penetrate 
Ene rin. | 

There are two more problems to,be coped with which do not have ap- 
parent solutions at this time. The first has to do with the superimposing 
of fringe information from different locations onto one scene beam line 
as illustrated in Figure 56. In case I the superimposing of the fringe 
change in Region A onto that in Region C and locating the fringe refer- 
ence line at pdnt C' is tolerable because the density field in Region A 
should be fairly constant and uniform. However in case ITI where the 
fringe change in all three regions are superimposed and located at point Gi. 
the answer is not readily apparent. If the plate and fin were of equal 
thicknesses then points A" and C'' would be the same and the error would 
be somewhat reduced. If, in addition, the thicknesses were made as thin 
as structurally possible, the error would be reduced to a minimum. Also 
it might be possible to integrate the model geometry into the computer pro- 


gram but this has not been attempted. 
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With the model at rotation angles between 0° and 90° scene beam de- 
fraction at the fin root and tip areas would also present problems as 
seen in Figure 57. A minimization of the problem could again be achieved 


by minimizing the fin and plate thicknesses. 


VI. CONCLUSIONS AND RECOMMENDATIONS 


The investigation, although not totally successful, has demonstrated 
the feasibility of using holographic interferometry to determine the 
flow field around a transparent model by looking through the model. The 
problems of model vibration and movement, of superimposing different 
fringe information on one beam, and of scene beam divergence through the 
fin root and tip regions must still be solved. The model vibrations and 
rotation to an angle of attack are attributable to the tunnel pressure 
fluctuations caused by tunnel leaks and to possible movement of the sting 
holder. It is felt that all three problems could be decreased to neglect- 
able effects by reducing the overall model size and by incorporating the 
model geometry into the computer pea 

The basic holographic arrangement was found to work quite well for 
this type of experiment. The holograms were generally high quality 
except when unfavorable tunnel flow conditions ance The use of 
circularly polarized light recommended by Okayama and Emori (14 | did not 
increase the hologram resolution appreciably but the method merits fur- 
ther considerations because of their excellent results. 

The data reduction process was found to be the rate-controlling step 
in the investigation due to the time and labor involved. Reducing the 
data from enlarged photographs of the interferogram saved some time and 


appeared to increase the data accuracy. Thesdata ceacter Of tae 3 fringe 
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was considered acceptable considering the fringe resolution in the holo- 
grams. It was felt that this could be reduced by using either a larger 
tunnel or a smaller model. With either of the changes it would be pos- 
sible either to use the free stream fringe pattern forward of the model 
as the reference conditions across the whole model or to use a vertical 
fringe pattern, since the leading edge Prandtl-Meyer expansion and fin 
shock would not block out the free stream fringes above the fin. It 
appears that the use of vertical fringes would also considerably reduce 
the data reduction time. 

With the use of Schlieren the flow network described by both Thomas 
N23, 24 | and Winkel mann [ 26, 27| was verified to exist. From the top 
view (90° model rotation) the fin shock and its fluctuations were ob- 
served and photographed through the flat plate plastic center section. 

The computer program, HOLOFER, was found to be quite capable of 
handling the type of flow field fringe data which would be generated in 
a complete analysis of this type. As pointed out before it is believed 


that the program should be modified to incorporate the model geometry. 
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Wind Tunnel Test Section with Clear Plastic Side Walls 


Figure l(a). 
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Figure 1(b). Wind Tunnel Test Section with Aluminum Side Walls 
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Figure 3. Holographic Arrangement Including the Laser Cooling and 
Aligning Equipment 
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Figure 4. Ruby Laser Power Supply and Cooling Unit 
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Figure 5. Ilologram Plate Holder and Mirrors on the Reverse Side 
of the Wind Tunnel 
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Hologram Platform Box Cover for Daylight Photography 


Figure 6. 
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Figure 7(a). Initial Fin-Flat Plate Model Used in the Experiment 
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Figure 7(b). Second Fin-Flat Plate Model Used in the Experiment 
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Figure 7(c). Details of the Second Fin-Flat Plate Model 
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Figure 7(d). Model Mounting in the Wind Tunnel Test Section 
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FIGURE §. CO-ORDINATE SYSTEM USED FOR THE INVERSION 
OF FRINGE NUMBER TO DENSITY 
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Figure 10. Desired Holographic Views in the Direct Fin-Root Method 
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Between G and 90° in the Direct Fin-Poot Mett.od 


50 





FROSTED 
PLEXIGLAS | 


Figure 12. 


STAINLESS STEEL 
REFLECTOR 







HOLOGRAM 
PLATE 


MIRROR OUTSIDE 
TUNNEL 


Schematic of the Model Used to Ottain Holographic Vicws 
Between 90° and 180° in the Direct Fin-Root Method 
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Figure 13. Holographic Viewing Angles Required in the Total Model 
Flow Method 
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Figure 14. Schematic of the Technique Used to Obtain Holograms in 
the Total tModel Flow !fethod 
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FIGURE 16. SPATIAL FILTERING TECHNIQUE FOR SELECTING 
PHOTOGRAPH OF CONSTANT ANGLE LINES OF 
LIGHT 
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FIGURE I7. LENSLESS PHOTOGRAPHIC TECHNIQUE USING A 
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Figure 18. Data Reduction Planes Desired to Describe the Density 
Field Down the Fin 
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Figure 20. Schematic of the Diffused Plastic Portion of the Model 
Required in the Direct Fin Flow Method 
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Figure 21. Interferogram Obtained in the Direct Fin Flow Method with 
the Model at O Rotation, Mach 2.84 and no Translation 
of Mirror, Me 
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Figure 22. Interferogram Obtained in the Direct Fin Flow Method with 
the Model at 0° Rotation, Mach 2.84 and a .006 inches 
Translation of Mirror, Mes Parallel to the Test Section 
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Figure 23. Interferogram Obtained in the Total Flow Method with the 
Clear Plexiglas Fin Model at O° Rotation, Mach 2.84 and 
a .0015 inches Horizontal Translation of the Diffuser 


Plate Parallel to the Test Section 
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Interferogram Obtained in the Total Flow Method with the 
Clear Plexiglas Fin Model at 0° Rotation, Mach 2.84 and 
a .0045 Vertical Translation of the Diffuser Plate 
Parallel to the Test Section 


Figure 24. 
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Figure 25(a). Lllustration of the Shock Network Around the Fin at 
Mach 2.84 with a 90° Model Rotation Angle Using a 


Horizontal Knife Edge Schlieren System 
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lllustration of the Shoct: Network Around the Fin at 
Mach 2.84 with a 90° Model Rotation Angle Using a 
Vertical Knife Edge Schlieren System 


Figure 25(b). 
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Figure 26. Illustration of the Model Flow Network at Mach 2.84, 0 
Model Rotation and a Plate Vibration of +.05 Angle of 
Attack Using a Schlieren System 
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Figure 27. ILIllustration of the Model Flow Network at Mach 2.84, 0° 
Model Rotation and No Plate Vibration Using a Schlieren 


system 
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Figure 40. Comparison of Fringe Data Obtained from Drawings of 
Interferogram Photographs 1 and 2 Aligned at Y' = 0.847 
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Figure 41. Comparison of Fringe Data Obtained from Photographic 
Enlargements of Interferogram Photographs 1] and 2 Aligned 
acy c= 0.847 
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Figure 42. Comparison of Fringe Data Obtained from Drawings of 
Interferogram Photographs 4 and 5 Aligned at Y' = 0.055 
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Figure 43. Comparison of Fringe Data Obtained from Photographic 
Enlargements of Interferogram Photographs 4 and 5 Aligned 
ate ) = 0.055 
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Figure 44. Comparison of Fringe Data Obtained from the Drawing and 
Photographic Enlargement of Interferogram Photograph 1 
Aligned at Y' = 0.847 
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Figure 45. Comparison of Fringe Data Obtained from the Drawing ard 
Photopraphic Enlargement of Interferopgram Photograpt 2 
Aligned at Y' = 0.647 
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Figure 46. Comparison of Fringe Data Obtained from the Drawing and 
Photographic Enlargement of Interferogram Photograph 3 
Aligned at Y’' » 0.451 


85 





1.4 


FROM PHOTO 4 
cp DRAWING 
1.2 se ENLARGEMENT 


ad 


0.8 


0.4 + 


0.2 





FRINGE NUMBER 


Figure 47. Comparison of Fringe Data Obtained from the Drawing and 
Photographic Enlargement of Interferogram Photograph 4 
Aligned ar Y' = 0.055 
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Figure 48. Comparison of Fringe Data Obtained from the Drawinp and 
Photographic Enlargement of Interferogram Photograph 5 
Aligned at Y' = 0.055 
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Figure 49. Comparison of Fringe Data Obtained from the Drawings and 
Photographic Enlargements of Interferogram Photoprapnrs l 
and 2 Aligned at Y' = 0.&47 
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Figure 50. Comparison of Fringe Data Obtained from the Drawings and 
Photographic Enlargements of Interferogram Photographs 4 
and 5 Aligned at Y' = (0.055 
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Fipure 51. Fringe Number Across the Fin in the Z = 0.387 Plane for 3 2 
Mach 2.84 as Determined from the Drawings of the Intericrorrem 
Photographs 
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Figure 52. Fringe Number Across the Fin in the Z = 0.387 Plane for 
§ = 0°, Mach 2.84 as Determined from the Photepraphic 
Enlargements of the Interferogram Photographs 
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Figure 53. Comparison of the Fringe Numbers Across the Fin in the 


2 = 0.387 Plane for = O°, Mach 2.84 as Determined 


from the Drawings and Photograrhic Enlargements of the 
Interferegram Photographs 
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Comparison of the Density Distributions Calculated by 
HOLOVLR for an Axisymmetric Case 
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Figure 55. Schematic of the Model Center Section Rotated to Illustrate 
the Loss of Fringe Information Due to Model Shadows on 
the Hologram 
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Figure 56. Schematic Illustrating the Problem of Different Fringe 


Information Being Superimposed on One Beam Caused by the 
Model Plastic 
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Figure 52. Schematic Illustrating the Scene Beam Refraction in the 
Fin Root and Tip Areas 
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APPENDIX A 
REDUCTION OF AN INTERFEROGRAM TO OBTAIN FRINGE SHIFT DATA 


The fringe shift reduction process was accomplished using two tech- 
niques. The first involved projecting the interferogram negative onto 
a sheet of white paper using a photo-enlarger. The light fringes offered 
the best contrast and were therefore traced out in Figures 30-34. In 
each drawing it was necessary to begin tracing the fringes above the fin 
and work towards the fin root since the transition across the fin tip 
determined the correct connection of the fringes across the fin leading 
edge shock. In order to determine the fringe change, one fringe line 
in the free stream region forward of the fin which appeared the straight- 
est and paralleled the majority of other fringe lines was selected. A 
Straight line, called the fringe reference line, was drawn over its 
centerline and extended to cross the y' axis. The remaining reference 
lines were then drawn parallel to the first and along the centerlines 
of the remaining free stream fringes. In reducing the drawings it was 
not realized until later that the free stream fringe patterns before and 
after the plate leading edge Prandtl-Meyer expansion differed consider- 
ably. This effect was taken into account later. 


' axis and 


A ruler scaled to 0.01 inches was then placed along the y 
the distances of the fringes and reference lines above and below the 
aligned y' plane were then recorded in Table I for Photograph 4. The fin 
width and length were then measured and the average values were recorded. 
The average fringe interval was determined by measuring the distance be- 


tween the first and last fringes used and dividing by the number of inter- 


vals. The fringe change was found by subtracting the fringe crossing 
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point. The fringe number was calculated by dividing the change by the 
average fringe interval. The reference line location was then normalized 
with respect to the measured fin height. Since the actual fin height 

was very close to one inch, the above number was considered to be the 
number of inches above or below the aligned plane. Thus the table com- 
puting the tunnel wall and grid refraction displacements in Appendix B 
could be entered to determine the actual reference line location with 
respect to the aligned plane of the drawing (see Figure 29). The loca- 


tions were then converted to the y' 


axis system by adding the normalized 
location of the aligned plane. After realizing that not all the refer- 
ence lines were referenced to the free stream density forward of the 
plate leading edge, the enlarged photographs (Figures 35-39) were checked 
against the drawings. It was found that the Prandtl-Meyer expansion 
caused a fringe num ber eranee of approximately 1.1. Consequently each 
reference line in the drawing was compared with the enlarged photograph 
and an appropriate percentage of the 1.1 fringe number was uSed as a 
correction (see Table I). The calculated fringe numbers had the cor- 
rection fringe number added to them while the y' locations were corrected 
by 


(1.1 Fringe no.) x (Fringe Interval ) 


y'corr = y'orig - (Fin Height) (A-1) 


The second reduction technique was to use enlarged photographs made 
from the interferogram negatives to obtain the fringe change. The fringe 
lines were first traced over lightly with a pencil and then verified 
against the other photographs to ensure correct tracing. A datum fringe 
reference line was chosen as before and drawn. The remaining reference 
lines for the upper free stream fringe lines forward of the Prandtl-Meyer 


expansion were then drawn parallel to the datum. When it became impossible 
she, 





to use fringe lines forward of the expansion, the reference lines were 
then drawn along the centerline of the fringes between the expansion 

and the fin. The same ruler was used to obtain the fringe crossing 
points, the reference crossing points, fin measurements, and the average 
fringe interval and they were recorded in Table 2. The reference line 
crossing points were then corrected to free stream conditions, using a 1.1 
Fringe number correction for those referenced to the pattern between the 
expansion and fin leading edge. The fringe Annee and reference fin 
locations were aleulatee as before. This method was considered more 
accurate because the reference lines and fringe lines could easily be 
rechecked for accuracy and corrected in the event that a fringe line 

was traced incorrectly or a reference line was misaligned. The accuracy 
was directly proportional to._the hologram resolution which was not true 
for the drawings since the reference lines are drawn parallel to handc- 


drawn fringe lines. 
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APPENDIX B 
CALCULATION OF TUNNEL WALL AND GRID PLASTIC REFRACTION CORRECTION 


In the interferogram photographs used to obtain the fringe change, 
every point off the alignment axis will be slightly distorted due to the 
plastic tunnel wall and grid. This effect is illustrated in Figure 19. 
From Snells Law of Refraction, the angles of incidence and refraction are 
related by 


sin& =n sinB (B-1) 


where n is the index of refraction between plastic and air. Then can be 


Written 


6 = ee ( sine ) 
= sin io ae (B-2) 


Since the tunnel wall and grid have the same index of refraction, they 
can be considered on material with thickness t = ab in Figure 19. Then 
consider the height bd which can be written 

bd = t tans (B-3) 


The beam displacement, Ay, is then 


AY = bd - be = t tan& - t tanp (B-4) 
but tan & is 
tan A= Yopserved (B-5) 
L 


Combining Equation (B-2), (B-4), and (B-5) the beam displacement becomes 


AY = if observed > tan [ sin”! sing (B-6) 


The true location of the observed point is then 


= - B-/ 
erie v observed Ay ( ) 
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A FORTRAN computer program was written to generate a table giving 
the true locations versus the observed locations. In the program the 


constants and variables from the above equations were defined as 


Yobserved es 

Ye = YTRUE 
AY = DY 
A = ALFA 

L = L 

n = 


Since the computer cannot calculate Equations (B-2) and (B-6) as written, 
they were constructed by parts using such letters as AA, AB, etc. The 


program and tables are included in the next few pages. 
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APPENDIX C 


APPLICATION OF COMPUTER PROGRAM "'HOLOFER" 


The computer program is an adptation of the inversion first proposed 
by C. D. Maldonado | 9, TO": 11| and is designed to invert fringe numbers 
across a field to the density field. It can be operated in three dif- 
ferent modes as described below: 

(a) Mode 1 

Mode 1 is utilized as a self-test of the computer program. [ft 
can either generate its own input density field using Subroutine FUNCT 
or read in a density field through Subroutine FREAD. The program then 
generates the fringe array and inverts the array back to the original 
density field. This mode was utilized in the present investigation to 
determine the value of the scale factor,c& , required to obtain the 
correct density across the fin. 

(b) Mode 2 

This mode reads in irregularly spaced fringe data and generates 
the fringe array at regular intervals across the field using Subroutine 
SHEET. By specifying NCODE = 1, the fringe array can be generated by one 
of the functions in Subroutine FUNCT. Mode 2 was not utilized. 

(c) Mode 3 

Mode 3 reads in the fringe data at regularly spaced intervals and 
inverts the array to density data across the field. The Subroutine GARRAY 
calls Subroutine READ to read in the fringe data. The first two cards 
preceeding the fringe data provide the program with the fringe field size, 


location, and symmetry. 


5 





The following parameters were used in considering the symmetric field 


case: 
PARAMETER INPUT 
NOF Run Number 
IMAX 201 
JMAX 1 
IS YM 101 
JSYM 1 
IMS Z 
JMS 100 
Z O23a7 
XO 0.0 
YO O20 
PHISYM | 0.0 


References [ 3 | and [12] contain further details and applications 
of the computer program. A print-out of the program is included in the 


and 


next few pages of this appendix. 
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